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ABSTRACT 

Variability amplitudes larger than 1 magnitude over time-scales of a few tens of minutes 
have recently been reported in the optical light-curves of several blazars. In order to inde- 
pendently verify the real occurrence of such extremely violent events, we undertook an ob- 
servational study of a selected sample of three blazars: PKS 0048-097, PKS 0754+100, and 
PKS 1510—089. Possible systematic error sources during data acquisition and reduction were 
carefully evaluated. We indeed found flux variability at intra-night time-scales in all three 
sources, although no extremely violent behaviour, as reported by other authors, was detected. 
We show that an incorrect choice of the stars used for differential photometry will, under fairly 
normal conditions, lead to spurious variability with large amplitudes on short time-scales. 
Wrong results of this kind can be avoided with the use of simple error-control techniques. 
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1 INTRODUCTION 

One distinctive characteristic of active galactic nuclei (AGN) is 
the fact that their energy flux varies along the whole electromag- 
netic spectrum, spanning a wide range of time-scales. Indeed, 
flux variability is an often used criterion for AGN detection (e.g.: 
MushotzkvlEoO^ . Within the widely accepted canonical model, 
i.e., super-massive black-hole (SMBH) + accretion-disk (AD) + 
dusty torus (DT) + relativistic jets, different regions of the active 
nucleus are thought to contribute to the power emitted at different 
frequencies along the spectral energy distribution (SED). On the 
other hand, variability time-scales strongly constrain the sizes of 
the emitting regions, through light-travel arguments. Since, except 
for the Mpc-scale radio-jets and lobes, the other components re- 
main spatially unresolved to astronomical observations in the vast 
majority of AGNs, flux variability studies are thus a powerful tool 
to probe those innermost regions. 

While most AGNs are variable on time-scales of a few years, 
blazars, i.e., the subclass comprising BLLac objects and flat- 
spectrum radio-quasars (FSRQ), display both the largest ampli- 
tudes and shorter time-scales. At optical wavelengths, long-term 
variations, with amplitudes of ~ 2 to ~ 5 mag along a few 
years have been found through extensive monitoring in several 
objects, such as GC 0109+224 jCiprini et al]|2003h . AO 0235+64 
iRaiteri et alJl2005l) . and OJ287 <Oian & Taolbooa) . Optical fluc- 
tuations on time-scales spanning several days are also usually ob- 



* Based on observations collected at the Complejo Astronomico El Leon- 
cito, operated under agreement among CONICET, UNLP, UNC, and UNSJ, 
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served, although with smaller amplitudes (e.g., ~ 1 mag in 20 days 
for PKS 200 5-489. iDominici et alj|2004 ~ 1.7 mag in 10 days 
for 3C 454.3. lFuhrmann et alfcOOffl . 

Although photometric changes for blazars on very short time- 
scales were reported more th an 30 y ears ago (Am ~ 0.3 mag 
in less than 24 hs for BLLac, |RacineJ|l970), this phenomenon re- 
mained unrecognized by most astronomers until the advent of CCD 
detectors, when the existence of the now called intra-night opti- 
cal var iability (INOV) or microvariability was fi rmly established 
jMiller. Carini & Goodrichll989tlCarini. Miller & Goodrichll99rJ: 
ICarini et alj Il99ll Il992h . Soon it became clear that several 
sources could experience remarkably large intra-night fluctua- 
tions, amounting to several tenths of a magnitude in a few hours. 
This i s the case, among others, for PKS 0537— 441 tfa nzjetalj 
Il986t iHeidt & Wagne3ll996t iRomero. Cellone & Combill2 000b). 
3C371 lCarinietalJl998). and AO 0235+64 l Rabbett e et al J 19961 



Heidt & Waaner 1996; Noble & Miller 1996). For this last object, 



Romero, Cellone & Combi 1 2000a) found changes up to 0.5 mag 
within one night and ~ 1.2 mag between consecutive nights, 
through well-sampled V and R band light-curves. This was one of 
the most violent variability events ever observed at optical wave- 
lengths in any blazar. The statistical incidence of microvariabil- 
ity in different classes of AGNs has been studied by several au- 
thors iRomero. Cellone & Combll99lllGopal-Krishna et all2003t 
ISaear et alj2004IStalin et alJ2005l) . 

The observed optical flux in blazars originates in a part 
of the AD and the inner (pc-scale) portions of the jets. Disre- 
garding gravitational micro-lensing effects, which probably ap- 
ply t o a small subset of particular sources (e.g.: iNilsson et alJ 
1996), three broad classes of intrinsic models may explain op- 
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tical microvariability: a) hot-s pot mode ls, involving instabili- 
ties in the AD (Manealam & Wiita|[l993l) . b) shock-in-jet mod- 
els, based on a relativistic shock-fr ont interacting w i th in- 
homogeneities or bends in the jet jM arscher & GeaJ Il985t 
iRomerd 1 1995k iRomero. Combi & Vucetichl Il995l) . and c) geo- 
metrical models where the jet changes its orientation relative to 
the observer, thus changing its Doppler factor jMarscheJfl995 
iGopal-Krishna & Wiitalll993) . Violent phenomena as those above 
reported, probably rule out models based on AD instabilities, need- 
ing, instead, relativistic beaming with high Doppler factors. 

In recent years, several papers claimed the repeated detection 
of extremely violent variability events in some b lazars, with am- 
plitudes Am > 1 mag in a few tens of minutes |Bai et alJll998l : 
IXie et aljll999ll200lL l2002Jbl 120041 : iDai et aljEoOlh . For exam- 
ple, these authors reported a 2 mag variation in ~ 40 min for the 
highly polarized quasar (HPQ) PKS 1510-089; they also showed 
several sudden (At < 15 min) "dips" of ~ 0.9 mag in the light- 
curve of the EGRET blazar OJ 248 (0827+243). If confirmed, these 
extremely violent phenomena would require a complete reassess- 
ment of the mechanisms which are thought to be responsible of the 
energy generation in AGNs. As an illustration, let us mention that 
optical variability time-scales of a few tens of minutes would imply 
emitting regions smaller than the Schwarzschild radius for certain 
objects. 

I ' contradiction to IXie et"all s claims. IRomero et alj (2002) 
found no such extremely violent variability events in a sample of 
20 EGRET blazars observed with ~ 20 min time-resolution along 
two or more nights each, showing that the discrepancies between 
both works probably had their origin in differe nt methods for er- 
ror control. In particular. IRomero et alj i200 2) suggested that an 
inappropriate choice of the comparison and control stars used for 
differential photometry could result in spurious fluctuations in the 
differential light-curve. 

Since it is desirable to firmly establish which is the real mini- 
mum time-scale of blazar microvariability, we have undertaken an 
intensive monitoring campaign focused on three particular objects 
with previous claims of extremely violent flux variations, in or- 
der to check this behaviour independently. We have made a careful 
analysis of the error sources involved in the differential photometry 
following well-established procedures, thus producing light-curves 
in which the significance of any variation is quantitatively evalu- 
ated. This procedure allowed us to set on solid bases the microvari- 
ability behaviour of blazars, showing that most claimed extremely 
violent events, if not all, are very likely to be spurious results pro- 
duced by an inappropriate error handling. In Sect.|2|we give details 
on the selected objects and the observations, while in Sect. [3] we 
describe our methodology, with special attention to the statistical 
error analysis. We show our results in Sect. [4] Spurious variability 
results are exemplified in Sect. [5] We close in Sect. [6] with some 
recommendations to future observers. 



2 SAMPLE AND OBSERVATIONS 

In order to independently check the reported events of extremely vi- 
olent microvariability in blazars, we have selected three of the most 
varia ble object s acco rding to IXie et alj 's papers lXieetalJ ll999. 
1200 ll 12002 jbl l2004 IDai et all200 ll) .Their names, equatorial co- 
ordinates, redshifts, catalogued visual magnitudes, and classifica- 
tion are listed in columns 1 to 6 of TableQ respectively. Here fol- 
lows a short description of each blazar. 



PKS 0048— 097: The optical spectrum of this BLLac shows very 
faint e mission lines, making it s redshift determination rather un- 
certain iRector & Stockel200ll) . A historical (~ 30 years) B-band 
light-curve shows a 1.8 mag variation; in the visual band, flux 
changes up to A V — 2.7 mag have been recorded. A larger vari- 
ation is reported at infrared wavelengths, amounting to ~ 6 mag 

1 Fan& Linll999h . On time-scales of several months! lFalomo et alj 
19931) reported a AV ~ 0.9 mag variation. 
IXie et alj l20 02fa) reported variations up to AR — 0.32 mag in 30 
minutes (January 2001). 

PKS 07 54+100: This is anothe r BL Lac whose re dshift is still un- 
certain jFalomo & Ulrichl EoOO). Baumert ( 1980) reported varia- 
tions up to ~ 2 mag in its optical flux, over long time-scales. 
This object also displayed fast polarization variability, both at op- 
tical (P p Dt = 4 - 26%) and IR (P m = 4 - 19%) wave- 
ij|l980). A l/-band light-curve compiled 
^ 3 mag change in 10 years, with 
1 mag in about 1 year. 



lengths l Puschell & Steir 
by Fan & Lin 1 2000) shows a 
smaller v ariation s up to AV ~ 

1998) reported AB = 0.47 mag in 22 minutes, while 
20041) claimed AB — 0.56 mag and AR = 0.66 mag in 



Bai ct al 



Xie et al 



about 80 minutes. 

PKS 1510-089: This is a well-studied FSRQ , 
X-ray spectrum l Singh, Shrader & George! 1 1 99^ ' 



with a hard 
and a power- 



ful gamma-ray emission, de tected by EGRET fThompson et alj 
Il993t ISreekumar et aljfl996h . At radio frequen cies, it has shown 
fast, large amplitude flux changes lAller. Aller & Hodgd fl98lL 
lAller. Aller & Hughe j|l99d) . Significant optical variations were 
first reported by Lii 1 1972) over a ~ 5 years tim e-scale. Its his- 
torical light-curve since 1899 was reconstructed bv lLiller&LilleJ 
1 1975); it shows a long-term variation with a maximum range 
AB = 5.4 mag, including an outburst in 1948 , after wh i ch the 
source brightness faded by 2.2 mag in 9 days. iGhoshetail f 2000) 
report "irregular variability of this blazar on time-scales of days to 
weeks," with a AR ~ 0.5 mag brightening in 84 days. 
Extremely violent events, with the highest amplitudes and shortest 
time-scales, were repeatedly claimed for this object: AR = 0.65 
mag in 13 min <Xie et alJEoOlh . AR = 2.0 m ag in 42 mi n 
iDai et alj200lh . and AV = 1.68 mag in 60 min <Xie et alj20023) . 
A AR = 1.35 ma g "dip" lasting At = 89 minutes was reported 
bv lXie etalll2004l) in the light-curve of this blazar. 

These three objects were the targets of our monitoring cam- 
paign, using the 2.15 m "Jorge Sahade" telescope at CASLEO, Ar- 
gentina, equipped with a Roper-EEV 1340 x 1300 pixels CCD 
(gain: 2.3 electrons adu" 1 ; read-out noise: 7.6 electrons). A focal- 
reducer provided a ~ 9 arcmin diameter field, with a scale of 
0.67 arcsec pix -1 . The blazars PKS 0048-097 and PKS 15 10-089 
were followed during six consecutive nights in August 2004, while 
PKS 0754+100 was observed along 5 nights (with a one-night gap, 
due to bad weather) in January 2005. Atmospheric conditions were 
photometric for ~ 54% of our observations, with some cirrus 
and/or bad seeing during the remaining time. Moon illumination 
was always below 38% during the August 2004 run and below 69% 
during the January 2005 run. 

Each object was repeatedly observed alternating between an R 
(Cousins) and a V (Johnson) filter. Integration times ranged from 
150 to 300 s; thus, between 6 and 20 images of the same object and 
with the same filter were obtained each night. 

Science frames were bias-corrected and then flat-fielded using 
master bias and flat-field frames (one flat for each filter) obtained 
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Table 1. Selected objects. 



Name 




"2000.0 


^2000.0 


z 




Type 






(h min s) 


P ' ") 




(mag) 




PKS 0048- 


097 


00 50 41.3 


-09 29 05 


0.22 


17.4 


BL Lac 


PKS 07544 


-100 


07 57 06.6 


+09 56 35 


0.27 


14.5 


BL Lac 


PKS 1510- 


089 


15 12 50.3 


-09 06 00 


0.36 


16.5 


HPQ 



by averaging 25 individual frames. The IRAF 1 package ccdred was 
used for this purpose. 



3 DATA ANALYSIS 

3.1 Statistical error analysis 

The usual technique to obtain light-curves for AGNs (and many 
other astronomical objects) is differential photometry, i.e., the mag- 
nitude of the target is measured against that of a comparison star 
recorded on the same CCD frame. In this way, small fluctuations 
due to non-photometric conditions are cancelled, since all star-like 
objec ts on the field of view are equally affected tHowell & Jacobvl 
1986). A second star, measured against the same comparison star, 
is often used as a stability check. In this way, variability is eval- 
uated by comparing the dispersions of two light-curves: target- 
comparison star (or) on the one hand, and control star-comparison 
star (a) on the other. Assuming that the second light-curve should 
only be affected by instrumental variations (i.e., both stars are non- 
variable), a statistical cri terion is often used (e.g. Ijang & Milled 
ll997tlRomero et all 19991) by introducing a parameter C = ot/c 
and requiring C > 2.576 for the source to be considered as variable 
at a 99% confidence level. 

This technique is conceptually so simple that possible prob- 
lems arising from its misuse are sometimes looked-over, thus lead- 
ing to the feeling that differential photometry is almost immune to 
any error source. To prevent this, several works have made use- 
ful recommen dations t o obse rvers, based on firm statistical bases. 
For example, Howell 1 1989) showed that very small (optimum) 
apertures maximize S/N for point-source observations, when used 
along with CCD growth curves to obtain instrumental magnitudes. 
However, it was later shown that such small photometric apertures 
should not be used for the particular case of an AGN embedded 
in a bright host galaxy, since this can lead to spurious variabil- 
ity results when the seeing FWHM changes along the observations 
ICellone. Romero & CombifcOOd) . 

Regarding the correct c hoice of stars for comparison and con- 
trol, a key work is that by Howell, Warnock & Mitchell 1 1988, 
hereafter HWM88). These authors show that it is not sufficient to 
simply select non-variable field stars for that purpose; these stars 
should also closely match the target's magnitude (colour matching 
is shown not to be so important). If not, the measured dispersion of 
the target-comparison light-curve (or) will be different from that 
of the control-comparison light-curve (a), just from photon statis- 
tics and other random-noise terms (sky, read-out noise), even in the 
absence of any intrinsic variations in the target. 



1 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 



Since suitable stars are not a lways found (especially for high- 
Galactic latitude fields), HWM88 give detailed calculations to de- 
rive a corrective factor V which properly scales a in order to 
match the expected instrumental dispersion <r T (inst) of the target- 
comparison light-curve. The computation of T thus requires the 
knowledge of the relevant CCD parameters, as well as mean values 
of the sky brightness and magnitudes of target and stars. Following 
HWM88 (their eq. 13) this corrective factor can be written as 



r " - m 



/Vfj (N T + P)+ N%(N S1 + P) 



N§ 2 (N T + P) + N%(Ns2 + P) 



(1) 



where TV stands for total (sky-subtracted) counts within the aper- 
ture, while sub-indices T, SI, and S2 correspond to the target, com- 
parison star, and control star, respectively. The factor P takes into 
account common noise-terms, being P — n p i x (/V s k y + N^ ON ), 
where n p i x is the number of pixels within the aperture, N s k y is the 
sky level, and TVron is the read-out noise. Median values are used 
for objects and sky. 

Thus, using the scaled a, the confidence parameter is now 
re-written as the quotient between the observed target-comparison 
light-curve dispersion and its expected dispersion just from instru- 
mental and photometric errors: 

r = rv = 



°"T (INST) 



(2) 



In the ideal case when all three objects are of the same magnitude, 
then F = 1, and the original definition of C is recovered. 

Many AGN variability studies, although not explicitly ap- 
plying HWM88's method, use comparison and control stars with 
apparent magnitude s very close to th at of the target object (e.g., 
iRomero et all 1 9991: ISaear et alfeooj . thus ensuring T ~ 1. Their 
results are in this way trustworthy, since the variability confidence 
levels of all light-curves are properly estimated. However, it is a 
fact that, whenever extremely violent variations have been claimed, 
they resulted from differential photometry using comparison and 
control stars more than ~ 2 mag (and up to ~ 5 mag) brighter than 
the AGN, and without any dispersion scaling. This flawed proce- 
dure leads to a severe overestimation of the confidence parameter. 
In what follows, we will further illustrate this point with results 
from our own observations. 



3.2 Photometry 

We used the IRAF package apphot to obtain aperture photome- 
try for the three blazars and several isolated, non-saturated stars in 
each field. Aperture radii were set at 8 pix (5'.'4, i.e., between ~ 1.5 
and 2 times the seeing FWHM) in order to prevent against any un- 
wanted effect due to light from the host galaxy under varying see- 
ing conditions l Cello ne et all200fj|) . We then selected the two most 
suitable stars in each field to be used as comparison (SI) and con- 
trol (S2), by requiring them to be non- variable and with magnitudes 
as close as possible to that of the blazar. In fact, the best resu lts are 
obtained when S 1 is slightly brighter than the target <HWM 88). so. 
we tried to fulfil this condition, too. 

We show finding charts for the three fields in Fig. Q Note 
that we generally did not use known standard stars in the blazars' 
fields, since most of these are too bright. The except ions are star 
S2 in the field of PKS 0048-097, which is star #4 in lVillata et ail 
ll998t). and star S2 in the field of PKS 754+100, which is star D 
in lMiller. Mullikin & McGimsevI Jl983l) . 

After constructing the differential light-curves, they were 
checked for any suspicious data-points, such as sudden changes 
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Figure 1. Finding charts for the observed fields. Each one is ~ 9 arcmin across, with North up, and East, to the left. Our target blazars are marked with T, 
while Si are the stars used for differential photometry (see text). 



in the control light-curve, and/or in only one of the photometric 
bands. A few of these events, due to cosmic-ray hits, were found 
and corrected. In the next section, we discuss our results. 



4 RESULTS 

4.1 Inter-night results 

Results from the whole campaign are given in Table [2] where 
columns 1 to 6 are, respectively: object name, filter, dates of ob- 
servation, control-comparison light-curve dispersion (a), total time 
between the first and last data points (At), and scaled confidence 
parameter (C/F). In column 7 we state whether the blazar was 
variable or not during our observations, according to the adopted 
criterion. The last two columns give, respectively, the maximum 
variation amplitude (in magnitudes) along the whole campaign, 
and the total number of data points in each light-curve. All light- 
curves are graphically shown in Figs. [2j|4] Each figure shows the 
target-comparison light-curve in the upper panel, and the control- 
comparison light-curve in the lower panel against heliocentric Ju- 
lian Date. Note that the vertical axis scale is always the same, in 
both panels of all three figures. 

It is evident that the three blazars were variable, with very high 
statistical significances, at these inter-night time-scales. Maximum 
amplitudes reached about half a magnitude for PKS 15 10— 089, and 
somewhat smaller values for the other two objects. Note that there 
is a good agreement between variability parameters in the V and R 
bands. 

Individual error bars and control-comparison dispersions are 
larger for a few nights, when scattered moonlight and/or track- 
ing errors affected our photometry. However, this had no evident 
impact on our ability to assess the object's variability. On the 
other hand, the data corresponding to the fifth observing night for 
PKS 15 10— 089 are affected by a 0.065 mag zero-point shift, due to 
technical problems. We corrected the graph in Fig.|4]for this effect, 
but we did not consider those data for the inter-night analysis. 



4.2 Intra-night results 

Tables |5||4| and [5] summarise the intra-night variability results for 
our three targets. Columns 1 to 6 give date of observation, filter, 
control-comparison light-curve dispersion (a), time spanned by the 
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2004 Aug 08-13 
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Figure 2. Differential V band light-curve for PKS 0048-097, whole cam- 
paign. Upper panel: AGN —Si ; lower panel: Sa — Si, 



PKS 0754+100 
2005 Jan 14-19 



. 4 » 




Figure 3. Same as Fig.|5]for PKS 0754+100. 

observations (At), and scaled confidence parameter (C/F), respec- 
tively. Column 5 states whether the blazar was considered to be 
variable or not during each night. The last column gives the num- 
ber of data points. We also show, for each blazar, the V light-curve 
for the night when the largest variation was detected (Figs.|5J0. 
All three blazars displayed microvariability, with amplitudes 



© 0000 RAS, MNRAS 000, 000-000 



Extremely violent microvariability in blazars 5 

Table 2. Inter-night results 



Object 


Filter 


UT Date 


a 


At 


c/r 


Variable? 


A/71 ln ax 


Nr. 






(dd-dd/mm/yy) 


[mag] 


M 






[mag] 




PKS 0048 -097 


R 


08-13/08/04 


0.007 


122.50 


13.61 


yes 


0.31 


53 




V 




0.007 


122.50 


13.07 


yes 


0.29 


54 


PKS 0754+100 


R 


14-19/01/05 


0.008 


125.21 


10.18 


yes 


0.36 


74 




V 




0.007 


124.94 


12.97 


yes 


0.42 


74 


PKS 1510—089 


R 


08-13/08/04 


0.006 


122.24 


31.89 


yes 


0.61 


55 




V 




0.007 


122.28 


23.99 


yes 


0.57 


58 



0.4 - « PKS 1510-089 

^ 2004 Aug 08-13 

0.6 - ^ 

0.8 - 



o.8 :::;:::;:::;:::;:::;:: 
■J * * * V f % 

< 1.0 - 

335 336 337 338 339 330 331 

JD-2452000 

Figure 4. Same as Fig.|2|ior PKS 1510-089. 



Table 3. Intra-night results for PKS 0048-097 



UT Date 


Filter 


a 


At 


Variable? 


c/r 


Nr. 


(dd/mm/yy) 




[mag] 


[h] 








08/08/04 


R 


0.004 


1.44 


no 


1.35 


6 




V 


0.008 


1.44 


no 


0.76 


7 


09/08/04 


R 


0.005 


2.33 


no 


1.37 


8 




V 


0.007 


2.38 


no 


0.27 


8 


10/08/04 


R 


0.008 


5.11 


yes 


3.22 


10 




V 


0.009 


5.11 


yes 


2.89 


10 


1 1/08/04 


R 


0.007 


2.38 


no 


0.59 


8 




V 


0.005 


2.38 


no 


0.80 


9 


12/08/04 


R 


0.006 


1.66 


no 


2.56 


9 




V 


0.006 


1.66 


no 


1.71 


7 


13/08/04 


R 


0.004 


2.83 


yes 


4.82 


12 




V 


0.004 


2.78 


yes 


5.54 


13 



up to Am ~ 0.08 mag in 1 hour; however, each of them was clas- 
sified as non- variable for at least one night. This does not mean that, 
in such cases, the object's flux was completely constant; all we can 
say is that any possible variation was then below our confidence 
threshold. 

It is thus clear that no extremely violent behaviour was de- 
tected in any source along our whole campaign. The statistical sig- 
nificance of this result can be assessed as follows. Let us define the 
duty cycle (DC) for extremely violent microvariability as the frac- 
tion of the observing time for which the object displayed large am- 
plitude (Am > 0.5 mag), fast (At < 45 min) flux changes. From 
the works by Xie's group (see Sect. a duty cycle DC ~ 50% 



Table 4. Intra-night results for PKS 0754+100 



UT Date 


Filter 


<7 


At 


Variable? 


c/r 


Nr. 


(dd/mm/yy) 




[mag] 


[hs] 








14/01/05 


R 


0.004 


1.82 


yes 


9.66 


8 




V 


0.011 


1.27 


yes 


4.02 


8 


15/01/05 


R 


0.007 


3.91 


yes 


5.78 


10 




V 


0.006 


3.91 


yes 


7.23 


8 


16/01/05 


R 


0.007 


4.75 


yes 


4.24 


19 




V 


0.008 


5.02 


yes 


4.58 


20 


17/01/05 


R 


0.010 


5.33 


yes 


3.01 


20 




V 


0.005 


5.33 


yes 


6.11 


20 


19/01/05 


R 


0.007 


4.56 


no 


1.94 


17 




V 


0.006 


4.51 


no 


1.81 


18 



Table 5. Intra-night results for PKS 1510-089 



UT Date 


Filter 




At 


Variable? 


c/r 


Nr. 


(dd/mm/yy) 




[mag] 


[hs] 








08/08/04 


R 


0.002 


2.38 


yes 


12.88 


6 




V 


0.004 


2.45 


yes 


6.68 


6 


09/08/04 


R 


0.004 


3.74 


yes 


3.38 


12 




V 


0.005 


3.74 


yes 


3.15 


12 


10/08/04 


R 


0.006 


3.17 


yes 


3.91 


10 




V 


0.004 


3.17 


yes 


3.75 


10 


1 1/08/04 


R 


0.004 


2.38 


no 


2.16 


11 




V 


0.008 


3.77 


no 


1.01 


13 


12/08/04 


R 


0.013 


3.17 


no 


1.10 


6 




V 


0.008 


3.24 


no 


1.16 


6 


13/08/04 


R 


0.005 


3.38 


no 


2.42 


10 




V 


0.009 


3.43 


no 


1.62 


11 



is inferred; if we accept this number, the fact that we did not de- 
tect any such extremely violent event along 17 observing sessions 
would have a very low probability (~ 10~ ). 

This result strongly implies that extremely violent microvari- 
ability reports should be carefully analysed, disentangling real flux 
changes in the source from systematic errors. In the next section 
we show how spurious results can be obtained from wrong error 
handling. 



5 SPURIOUS VARIABILITY 

As an illustrative example of the effects of a bad choice of compar- 
ison and control stars, let us consider our data for PKS 1510—089 
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-0.15 
2 -1.15 



PKS 0048-097 
2004 Aug 13 




UT [h] 



2004 Aug 10 
Field star 
(Spurious variability) 




UT [h] 



Figure 5. Differential V band light-curve for PKS 0048-097, on the night 
of Aug 13, 2004. Upper panel: AGN -Si ; lower panel: Sa— Si. 



PKS 0754+100 
2005 Jan 15 



UT [h] 

Figure 6. Differential V band light-curve for PKS 0754+100, on the night 
of Jan 15, 2005. Upper panel: AGN —Si ; lower panel: Sa— Si. 



PKS 1510-089 
2004 Aug 08 



J s s 



UT [h] 

Figure 7. Differential V band light-curve for PKS 1510-089, on the night 
of Aug 8, 2004. Upper panel: AGN -Si ; lower panel: S 2 -Si . 



Figure 8. Differential V band light-curve for a faint non- variable star in the 
field of PKS 1510-089, on the night of Aug 10, 2004. Upper panel: S 5 - 
S6; lower panel: S7— Sg. The apparent extreme variability is just a spurious 
result. 



on the night of Aug. 13, 2004. Our originally selected stars approx- 
imately follow the prescriptions given in Sect. 13.11 AV(T — Si) = 
0.90, AV(Sa — Si) = 0.96; i.e., comparison star not much brighter 
than target AGN and control star slightly fainter than the AGN. 
With these stars, we obtain a confidence parameter C = 1.55, 
and a scaled confidence parameter ^ = 1.62; i.e., the blazar is 
classified non-variable in both cases. We now use two significantly 
brighter stars, shown as S3 and S4 in Fig.Q AV(T — S3) = 1.95, 



AVYS4 - S3) = 0.58 (S 3, S 4 



star #4, sta r #6 in lRaiteri et alJ 
119981 = star A, star B in lVillata et alJll997l respectively). This 
choice gives a confidence parameter C — 3.72, implying that the 
blazar should be considered as variable during that night. However, 
after correcting for the large flux difference between the objects, a 
scaled confidence parameter § = 1.32 is obtained, thus classifying 
the target as non- variable. 

The preceding analysis shows that, although a certain level of 
intrinsic variability may be present in the target, its significance is 
severely overestimated when stars ~ 2 mag brighter than the tar- 
get are chosen for comparison and control. We now want to test 
if dramatic, although spurious, variability events may be produced 
when still brighter stars are used. Since any star more than ~ 3 
mag brighter than our targets was always saturated on our images, 
we selected a faint, non-variable star in the field of PKS 15 10— 089 
to illustrate this effect. This star is labelled as S5 in Fig. Q and 
we used stars and S7 as compari son and control, respectively 
(S6 = star #5 in lRaiterietalJll998h . Hence, we reproduce a sit- 
uation where the target is ~ 5 mag fainter than the comparison 
star, while this last, in turn, is ~ 1.4 mag brighter than the control 
star. The result of this completely inappropriate choice of stars is 
shown in Fig. [8] in about one hour, the source seems to brighten 
by ~ 0.4 mag, then returning to its original flux level; later, a 
spectacular "outburst" brightens the object by ~ 1.2 mag in about 
35 minutes. Note the stability of the control - comparison light- 
curve; without any dispersion scaling, the variability confidence 
parameter is C — 24.0, thus giving the impression of a highly 
significant flux variation. However, the scaled confidence param- 
eter is just jr = 1-0, clearly establishing that no significant vari- 
ability is present in the data, with any fluctuation in the object - 
comparison light-curve being the result of photometric errors dif- 
ferently affecting stars of very disparate magnitudes. In this par- 
ticular case, both "outbursts" coincide with sudden changes in at- 
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mospheric transparency, due to passing cirrus. Several factors con- 
tribute to make errors larger during such events: a lower flux level 
from the star, higher sky level and its associated rms, changes in 
the seeing FWHM, etc. Brighter stars are relatively less affected by 
these effects than the faint target, thus leading to the apparent vari- 
ations in the differential light-curve. Under different observational 
conditions and with different photometric techniques, it is likely 
that spurious "dips" instead of "outbursts" may be produced in the 
light-curve. 

These results underscore the absolute necessity of using the 
method described in HWM88 whenever suitable comparison and 
control stars cannot be found. Not doing so will very probably lead 
to spectacular, although completely false, results. 

In some cases, ex tremely vi olent variability events have been 
claimed to be periodic IXie et all2 004). This is not surprising, since 
the observations are periodically repeated, and, if the same inade- 
quate photometric techniques are used, similar spurious results will 
be obtained. These systematic errors have generated some com- 
pletely ill-motivated theoretical models (Wu et al. 2005). 



6 SUMMARY AND CONCLUSIONS 

Repeated claims for the detection of extremely violent op tical vari- 
ability in blazars have been raised up in recent years IXie et all 
ll999ll200lLf2002jbLl2004lDai et all200ll) . These claimed events 
are characterised by fast flux changes Am > 0.5 mag in a 
few tens of minutes, and reaching up to, for example, a 2 mag 
variation in ~ 40 min reported for PKS 1510—089. However, 
other studies have found that the typical minimum time-scale for 
such large-amplitude vari ations in bl azars is of several hours or 
still larger (e.g., Ro mero et alj |2002j). We have thus undertaken 
an observational campaign, targeting the blazars PKS 0048—097, 
PKS 0754+ 100, and PKS 1510-089, devised to shed light on this 
controversy. This paper presents its results, showing that, although 
microvariability was clearly detected in our three targets, no ex- 
tremely violent optical variability event was detected along 110 
hours of observation. The largest fast flux variations we detected, 
instead, amount to < 0.1 mag in about one hour. 

We show that this discrepancy is most likely due to systematic 
errors introduced during the observations and photometry. In par- 
ticular, the use (without any correction) of stars much brighter than 
the target for differential photometry, directly leads to an overes- 
timation of the significance of any detected variability. Moreover, 
under certain specific conditions, it easily gives place to spurious 
variability closely reproducing extremely violent microvariability 
events as those reported in the papers by Xie and coworkers. 

The following recommendations should thus be followed in 
order to prevent against spurious variability results: 

The target object must be neither underexposed, nor saturated 
on all science frames. Differential light-curves should be made us- 
ing comparison and control stars as close in magnitude as possible 
to the target. Published standard stars in blazar fields are usually 
too bright for this purpose; they should be used just for calibration 
to the standard system, through a few short exposure-time frames. 

If no suitable stars can be selected for differential photometry, 
brighter (or fainter) stars may be used, provided that the variance of 
the control-comparison light-curve is prope rly scaled . This should 
be done following the meth od present ed in HWM88 (in fact, it is 
always recommended to use HWM88 method). 

Any remarkable flux change should be critically verified, 
looking for cosmic-ray hits, sudden changes in atmospheric trans- 



parency, or any instrumental effect (see also lCellone et all20 00 for 
spurious variability induced by seeing FWHM changes). 

We conclude by saying that a critical evaluation of past and fu- 
ture claims for extremely violent microvariability events in blazars 
is needed before any radical revision of blazar models be required. 



ACKNOWLEDGMENTS 

This work received financial support from ANPCyT (PICT 03- 
13291 BID 1728/OC-AR) and CONICET (PIP 5375). Additional 
support from IALP (UNLP-CONICET) is also acknowledged. We 
thank Cecilia Farina for help with the observations, and CASLEO 
staff members A. de Franceschi, R. Jakowczyk, and P. Ostrov for 
their skillful assistance at the telescope. We wish to thank the 
anonymous referee for his/her very constructive comments. 



REFERENCES 

Aller H. D., Aller M. F, Hodge P. E., 1981, AJ, 86, 325 
Aller M. F, Aller H. D., Hughes P. A., 1996, in Miller H. R., 
Webb J. R., Noble J. C, eds, ASP Conf. Ser. 110: Blazar Con- 
tinuum Variability Centimeter-wavelength Flux and Polarization 
Variability as a Probe of the Physical Conditions in AGN (I), 
p. 193 

Bai J. M., Xie G. Z., Li K. H, Zhang X., Liu W. W., 1998, A&AS, 
132, 83 

Baumert J. H, 1980, PASP, 92, 156 

Carini M. T, Miller H. R., Goodrich B. D., 1990, AI, 100, 347 
Carini M. T, Miller H. R., Noble J. C, Sadun A. C, 1991, AI, 
101, 1196 

Carini M. T, Miller H. R., Noble J. C, Goodrich B. D., 1992, AI, 
104, 15 

Carini M. T, Noble I. C, Miller H. R., 1998, AI, 1 16, 2667 
Cellone S. A., Romero G. E., Combi I. A., 2000, AI, 1 19, 1534 
Ciprini S., Tosti G, Raiteri C. M., Villata M., Ibrahimov M. A., 

Nucciarelli G, Lanteri L., 2003, A&A, 400, 487 
Dai B. Z., Xie G. Z., Li K. H, Zhou S. B., Liu W. W., liang Z. I., 

2001, AI, 122, 2901 
Dominici T. P., Abraham Z., Teixeira R., Benevides-Soares P., 

2004, AI, 128, 47 
Falomo R., Ulrich M.-H, 2000, A&A, 357, 91 
Falomo R., Bersanelli M., Bouchet P., Tanzi E. G, 1993, AI, 106, 

11 

Fan I. H., Lin R. G, 1999, ApIS, 121, 131 
Fan I. H., Lin R. G, 2000, Apl, 537, 101 
Fuhrmann L. et al., 2006, A&A, 445, LI 

Ghosh K. K, Ramsey B. D., Sadun A. C, Soundararajaperumal 

S., 2000, ApIS, 127, 11 
Gopal-Krishna, WiitaP. I., 1992, A&A, 259, 109 
Gopal-Krishna, Stalin C. S., Sagar R., Wiita P. I., 2003, Apl, 586, 

L25 

Heidt I., Wagner S. I., 1996, A&A, 305, 42 

Howell S. B., 1989, PASP, 101, 616 

Howell S. B., lacoby G. H, 1986, PASP, 98, 802 

Howell S. B., Warnock A. I., Mitchell K. I., 1988, AI, 95, 247 

lang M., Miller H. R„ 1997, AI, 114, 565 

LillerM. H., Liller W., 1975, Apl, 199, L133 

Lii P. K, 1972, AI, 77, 829 

Mangalam A. V, Wiita P. I., 1993, Apl, 406, 420 



© 0000 RAS, MNRAS 000, 000-000 



8 Cellone et al. 



Marscher A. P., 1992, in Zensus J. A., Wagner S. J., eds, Physics 
of Active Galactic Nuclei Emission Models for Compact Jets. 
Springer- Verlag, Heidelberg, p. 510 
Marscher A. P., Gear W. K., 1985, ApJ, 298, 1 14 
Miller H. R., Mullikin T. L., McGimsey B. Q., 1983, AJ, 88, 1301 
Miller H. R., Carini M. T., Goodrich B. D., 1989, Nature, 337, 
627 

Mushotzky R., 2004, How are AGN Found?. ASSL Vol. 308: Su- 
permassive Black Holes in the Distant Universe, p. 53 

Nilsson K., Charles P. A., Pursimo T., Takalo L. O., Sillanpaa A., 
Teerikorpi P., 1996, A&A, 314, 754 

Noble J. C, Miller H. R., 1996, in Miller H. R., Webb J. R., Noble 
J. C, eds, ASP Conf. Ser. 1 10: Blazar Continuum Variability Ex- 
treme Optical Microvariability in Three Selected Blazars. p. 30 

Puschell J. J., Stein W. A., 1980, ApJ, 237, 331 

Qian B., Tao J., 2003, PASP, 115, 490 

Rabbette M., McBreen S., Steel B., Smith N., 1996, A&A, 310, 1 
Racine R., 1970, ApJ, 159, L99 
Raiteri C. M. et al., 2005, A&A, 438, 39 

Raiteri C. M., Villata M., Lanteri L., Cavallone M., Sobrito G, 

1998, A&AS, 130, 495 
Rector T. A., Stocke J. T., 2001, AJ, 122, 565 
Romero G. E., 1995, Ap&SS, 234, 49 

Romero G. E., Combi J. A., Vucetich H., 1995, Ap&SS, 225, 183 
Romero G. E., Cellone S. A., Combi J. A., 1999, A&AS, 135, 477 
Romero G. E., Cellone S. A., Combi J. A., 2000a, A&A, 360, L47 
Romero G. E., Cellone S. A., Combi J. A., 2000b, AJ, 120, 1 192 
Romero G. E., Cellone S. A., Combi J. A., Andruchow I., 2002, 
A&A, 390, 431 

Sagar R., Stalin C. S., Gopal-Krishna G, Wiita P. J., 2004, MN- 
RAS, 348, 176 

Singh K. P., Shrader C. R., George I. M., 1997, ApJ, 491, 515 
Sreekumar P. et al., 1996, ApJ, 464, 628 

Stalin C. S., Gupta A. C, Gopal-Krishna, Wiita P. J., Sagar R., 

2005, MNRAS, 356, 607 
Tanzi E. G. et al., 1986, ApJ, 311, L13 
Thompson D. J. et al., 1993, ApJ, 415, L13 
Villata M. et al., 1997, A&AS, 121,119 

Villata M., Raiteri C. M., Lanteri L., Sobrito G., Cavallone M., 

1998, A&AS, 130, 305 
Wu J., Zhou X., Peng B., Ma J., Jiang Z„ Chen J., 2005, MNRAS, 

361, 155 

Xie G. Z., Li K. H., Zhang X., Bai J. M., Liu W. W., 1999, ApJ, 
522, 846 

Xie G. Z., Li K. H., Bai J. M., Dai B. Z., Liu W. W., Zhang X., 

Xing S. Y., 2001, ApJ, 548, 200 
Xie G. Z., Liang E. W., Zhou S. B., Li K. H., Dai B. Z., Ma L., 

2002a, MNRAS, 334, 459 
Xie G. Z., Zhou S. B., Dai B. Z., Liang E. W., Li K. H, Bai J. M., 

Xing S. Y., Liu W. W., 2002b, MNRAS, 329, 689 
Xie G. Z., Zhou S. B., Li K. H., Dai H„ Chen L. E., Ma L„ 2004, 

MNRAS, 348, 831 



